Critical Point of a Weakly Interacting Two-Dimensional Bose Gas 
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We study the Berezinskii-Kosterlitz-Thouless transition in a weakly interacting 2D quantum Bose 
gas using the concept of universality and numerical simulations of the classical |'!/;|'*-model on a lat- 
tice. The critical density and chemical potential are given by relations ric — {mT/2-Kh^) ln{^h^ /mU) 
and fj,c = (mrfZ/Trft^) In(^p7i^/ml7), where T is the temperature, m is the mass, and U is the ef- 
fective interaction. The dimensionless constant ^ = 380 ± 3 is very large and thus any quantitative 
analysis of the experimental data crucially depends on its value. For our result is — 13.2 ±0.4. 
We also report the study of the quasi-condensate correlations at the critical point. 



An accurate microscopic expression for the criti- 
cal temperature of the Berezinskii-Kosterlitz-Thouless 
(BKT) transition [|| has been a weak point of the the- 
ory of weakly interacting two-dimensional Bose gas. The 
theory of Ref . |^ (see also |^,^ and analysis below) , sug- 
gests that the critical density of the BKT transition in 
the weakly interacting system reads (we set h = 1) 



2tt mil 



(1) 



However, the value of ^ cannot be obtained within stan- 
dard analytical treatments since ^ is related to the system 
behavior in the fluctuation region where the perturbative 
expansion in powers of U does not work. With unknown 
^, one finds Eq. (^ rather inaccurate unless mil is expo- 
nentially small. Moreover, as we will find in this Letter, 
the value of ^ is very large: ^ « 380. This means that for 
all experimentally available up to date (quasi-)2D weakly 
interacting Bose gases the quantitative analysis of 
the data for the critical ratio ric/Tc requires a precise 
value of ^. In the system of spin-polarized atomic hydro- 
gen on helium film the value of mU is of order unity 
in the recently created quasi-2D system of sodium 
atoms B, mU is of order 10~^, according to the formula 
of Ref. 0. 

To quantitatively describe the BKT transition in a 
weakly interacting Bose gas, it is sufficient to solve 
a classical-field j-i/'l^-model with the effective long- wave 
Hamiltonian B 



(2) 



where fi' is the chemical potential, and tp is the classical 
complex field. 

In this Letter, we first discuss the origin of the relation 
(^ in the limit of small U, and how quantum and clas- 
sical models relate to each other. Then we present our 
numeric results (for the critical density, critical chemi- 
cal potential, and quasi-condensate correlations at the 
BKT point) obtained by simulating the critical behavior 



of the 2D I -01 ''-model on a lattice using recently devel- 
oped Worm algorithm |^ for classical statistical models. 
In particular, we show that the quasi-condensate correla- 
tions are very strong at Tc , in agreement with the experi- 
mental observation in the spin-polarized atomic hydrogen 
Q and quantum Monte Carlo simulations 

A simple dimensional analysis of the Hamiltonian (||) 
allows to write a generic formula for the critical point in 
a weakly interacting 2D j-^l^-model. The routine itself 
is completely analogous to that in the 3D case (see, e.g., 
JlO| , p| ), but final results naturally reflect the specifics of 
the 2D case. 

We begin with introducing the mode-coupling momen- 
tum, fee, that characterizes the onset of strong non-linear 
coupling between the long- wave harmonics of ^(r) (har- 
monics with k ^ kc are almost free). This momentum is 
just the inverse of the healing length, or vortex core ra- 
dius, Tc ||l|. We denote by h the contribution to the total 
density due to strongly coupled harmonics, and introduce 
the renormalized chemical potential 



(3) 



by subtracting the mean field contribution of non- 
interacting high-momentum harmonics. Here = 
( iV'kP), and (. . .) stands for the statistical averaging. 

An estimate for h follows from the Nelson-Kosterlitz 
formula 



2mT 
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since it is intuitively expected that n ~ n^. An indepen- 
dent estimate of the parameters of the fluctuation region 
is obtained by considering when all three terms in Eq. (|^) 
are of the same order: 



kl/m ^ ~ fiU 



(5) 



and relating h ~ X)fe<fc '^fc ~ ^c^fec to the renormal- 
ized chemical potential by using T/|/i| in place of the 



1 



occupation number n^^. By definition, kc separates 
strongly coupled and free harmonics, and thus Uk^ ~ 
T/[k1/2m — jl] ^ T/\il\. The final order-of-magnitude 
estimates read (at T = Tc) 



U 



fi ^ mT , 
fee ~ m{UTfl'^ 
jl ^ U mT . 



(6) 
(7) 
(8) 



We are now in a position to derive Eq. (|l|) for the crit- 
ical density. In 2D the main contribution to the integral 



n = / nkd?k/{2TTY 
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comes from large momenta between kc and some ultra- 
violet scale fc*. The value and physical meaning of fc, 
depend on the model. For classical lattice models fc, is 
given by the inverse lattice spacing; in the continuous 
quantum system fc* ~ V mT is the thermal momentum. 
Ax, kc ^ k <^ k^ we have nj. « 2mT jk^, and thus can 
write 



"^x^iCklikD 



(10) 



where C is some constant. Critical density, Eq. (|^), 
for the quantum Bose gas is obtained by substituting 
Ckljkl = imT/m^UT = ^/mU. 

The dependence on the ultraviolet cutoff is associated 
with the properties of ideal systems only, while the long- 
wave behavior of all weakly-interacting |^/;| ^-theories is 
universal. This fact allows one to relate results for dif- 
ferent models by adding and subtracting non-interacting 
contributions, i.e., up to higher order corrections in U 
the difference between models A and B is given by 



(n! 



(^) 
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) = !h 



(ideal A) 



^(idcal i3)]rf2^/(27r)2. 



In 



what follows, the reference system A will be the classical 
lattice model with lattice spacing a, and our results are 
analyzed using 



flat) _ 



mT 



In 



A 



27r m^a^UT 
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The actual system of interest is the quantum Bose gas, 
so we add and subtract the corresponding ideal-system 
contributions to get 



In 



A 



1 



S^ma-^T 2TrmT 



BZ 



T d^k 



d^k 
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where BZ means that the first integral is over the Bril- 
louin zone, and £'(k) is the dispersion law for the ideal 
lattice model such that E{^^ ^ 0) ^ k'^ /2m. [The diver- 
gences of the two integrals in Eq. (|l^) at fc ^ compen- 
sate each other.] 

Our simulations were done for the simple square lattice 
Hamiltonian 



H=Y. [i5(k)-^]|^,|2 + _^|^,|4, (13) 

where tpk is the Fourier transform of the complex lattice 
field ^Ji^ and 

E{\i) = {l/ma^)[2 - cos(fc^a) - cos(fcj^a)] (14) 

is the tight-binding dispersion law. With this dispersion 
relation the r.h.s. in (|l2|) can be evaluated analytically 
and we obtain the "conversion" formula 



^ = A/16 



(15) 



Since final results for dimensionless constants do not de- 
pend on TO, T, and a, in numerical simulations we set 
a = 1, r = 1, and m — 1/2 for convenience. 

The above consideration for the critical density can 
be readily generalized to the critical chemical potential, 
with the result 



mTU 
In— 

TT mU 



(16) 



First, we notice that Eq. ( p^ ) immediately follows form 
Eqs. (||) and Eq. (|^) because the mean- field term is pro- 
portional to —{mUT/T:)ln(mU) (we actually deal with 
exactly the same integral). Since the renormalized value 
jl is universal, to account for the difference between the 
classical and quantum models one has to add and sub- 
tract mean-field contributions dominated by the ideal 
behavior. Thus, if the classical model is analyzed us- 
ing /ic = {mTU/TT)ln[A^/m'^a^UT], one has to apply 
Cm = Af^/16 to get the quantum result, Eq. (^6|). 

We now turn to our numerical procedure. To simulate 
the grand-canonical Gibbs distribution corresponding to 
the Hamiltonian (|l3|), we employ the Worm algorithm 
(see Ref. for the description) that has demonstrated 
its efficiency for the analogous problem in 3D ||ll|. The 
formal criterion of the critical point for the finite-size 
system is based on the exact (Nelson-Kosterlitz) rela- 
tion (^): We say that the system of linear size L is at 
the critical point, if its superfluid density, ns{L), satis- 
fies ns{L) — 2mT/n. [The superfluid density has a di- 
rect estimator in the Worm algorithm via the statistics 
of winding numbers [^ , and its autocorrelation time does 
not suffer from critical slowing down.] 

The finite size scaling of nc{L) is well known from the 
Kosterlitz-Thouless renormalization group theory Q] 



nc{L) 



A' mT 



In^ [A" Lm{UTy 



/2] 
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where A' and A" are dimensionless constants. A sim- 
ilar relation applies also to the critical chemical poten- 
tial. Equation (^7|) was used for the finite-size scaling 
analysis. We found that instead of extrapolating data 



2 



for each value of U to the L ^ oo limit independently, 
a much more efficient procedure is to perform a joint 
finite-L and finite-?/ analysis. To this end we heuristi- 
cally introduce parameters accounting for non-universal 
finite-C/ corrections by adding linear in U terms to each 
of the three of the dimensionless constants: A — > A+BU, 
A' ^ A' + B'U, and A" A" + B"U. We thus have 
six fitting parameters to describe all our data points 
p3[ . The data for n(.{U, L) and ^ic{U, L) are presented in 
Fig. 1. The fitting procedure yields A = (6.07±0.05)•10^ 
A^ = (211±6), which, according to Eq. (|l^), means that 

(18) 



C = 380±3, ^^ = 13.2 ±0.4 



The fit is extremely good — 20 points for the critical den- 
sity at C/ < 2.5 and Lm{UTy^^ > 15, each calculated 
with relative accuracy of order 10~^, are described with 
the confidence level of 62 %. 



nc(U,U 
1.1- 

1- 


▼ 
♦ 




T 


T - 0.01 

A - 0.05 

♦ - 0.125 
o - 0.25 
■ - 0.5 

H - 1. 

• - 2.5 

» - 5. 


0.9- 

0.8- 


o 

■ 






o 


0.7- 


• 

• 

O 


1 




1 





0.5 




ln-2[Lm(UT)"2 


U 

1.8 








» - 0.01 
• - 0.05 


1.6 










B - 0.125 
■ - 0.25 
» - 0.5 


1.4 








' - 1. 
' - 2.5 


1.2 








. - 5. 


1 










0.8 




1 




1 





0.5 




In"' [Lm(UT)'"] 



FIG. 1. Critical density and chemical potential for various 
coupling parameters and system sizes. Typical error bars are 
much smaller than symbol sizes. The dotted line is the fitting 
function described in the text. 

Experiments on helium films often report that the ratio 



ns{Tc)/ns{0) = ns{Tc)/nc = ^mTc/irnc is close to 0.75 
Jl^ , p^ . Our simulation predicts that this ratio is given 

by 



ns{Tc)/nc 



5.94 - ln(mt/) 



(19) 



and mil w 1.8 is required to describe helium films, pro- 
vided the small- [/ approximation may be pushed that far 
]l6t . We are not aware of the published data on the crit- 
ical chemical potential. [For helium and hydrogen films 
on substrates one has to shift /ic by the value of the ab- 
sorption energy (for the delocalized atom, in the case of 



helium film), fic 



— En 



{mTU/Tr)\n{^f,/mU). In 



thermal equilibrium this quantity can be readily mea- 
sured through the chemical potential of the bulk vapor.] 
In the absence of long-range order parameter, 2D sys- 
tems below Tc are characterized by the local correlation 
properties of the quasicondensate density, identical to 
those of a system with genuine condensate These 
properties refiect the specific structure of the V'-held: 



V'(r) 
^o(r) 



„i$(r) 



noe 



(20) 
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where the quasicondensate density hq may be considered 
as a constant, and ^pi is the Gaussian field independent 
of ipQ. Both experiment and model Monte Carlo sim- 
ulations [|j indicate that in 2D systems with mil ^ 1 
the quasicondensate correlations appear well above Tc 
and are pronounced at Tc- Below we show that this is a 
generic feature of weakly interacting ji/ij^-models. 

It is convenient to characterize the quasi-condensate 
properties by the correlator 



(22) 



The Gaussian component of the field obeys the Wick's 
theorem and does not contribute to Eq. (^2|). If, for a 
moment, by ipi we understand short-wave harmonics of 
"0, we conclude that only long- wave and strongly non- 
linear harmonics with the momenta k contribute to 
the correlator Q, i.e. Q ~ n^. Thus, we expect that all 
weakly interacting |'0|'*-models satisfy 



Q 



(T = Tc) (23) 
is a universal constant. 



in the limit of small U, where C* 
By definition, hq — \fQ. 

The finite-size and small-J7 analysis of the data for 
Q{U, L) was done in complete analogy with previously 
discussed cases of nc{U,L) and fic{U,L) (see Ref. |^). 
We found that 



a = 1.30 ±0.02 



(24) 



The ratio between no(r — Tc) and Uc describes how pro- 
nounced are the quasicondensate correlations in the Bose 
gas at the BKT point: 



3 



7.16 



27r\/ C* 



\n{(/mU) 5.94 + ln(l/TOf7) 



(25) 



We see, that it is of order unity unless mil is exponen- 
tially small. Another interesting ratio is 



no TT\/CU 
n. 2 



1.79 



(26) 



which is interaction independent and shows that the su- 
perfluid density is substantially smaller than the quasi- 
condensate density at T^. 



m(UT)'' 



FIG. 2. Quasicondensate correlations as a function of sys- 
tem size. The dotted line is to guide the eye. 

Finally, we would like to derive an accurate estimate 
for the mode-coupling radius re- In an ideal system 
(3 = 0. Hence, Q{L) should decrease with decreasing L, 
and for system sizes L ~ Tc it has to drop significantly 
from its thermodynamic value. We rather formally de- 
fine Tc from Q{L = Vc) ~ Q{L oo)/2, and from Fig. 2 
obtain 



Tc « 2/TO([/r)i/2 



(27) 



We conclude by noting that Nelson-Kosterlitz formula 
(^ and Eqs. (|l]), (p^), and ( p3| ) constitute a complete 
set of equations which allow to fully determine system 
parameters from measurements with independent cross- 
checks. We are not aware of another study were dimen- 
sionless constants ^, and C* were determined with 
high precision. 
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